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\Q ' ABSTRACT 

<n : 

■ 

j> ■ We discuss the effects of synchrotron self-absorption in the prompt emission 

from the reverse shock of GRB afterglows in a dense environment, such as the 
wind of a stellar progenitor or a dense ISM in early galaxies. We point out that, 
when synchrotron losses dominate over inverse Compton losses, the higher self- 
absorption frequency in a dense environment implies a bump in the reverse shock 
emission spectrum, which can result in a more complex optical/IR light curve 
than previously thought. This bump is prominent especially if the burst ejecta is 

q . highly magnetized. In the opposite case of low magnetization, inverse Compton 

+3 . losses lead to a prompt X-ray flare. These effects give a possible new diagnostic 

for the magnetic energy density in the fireball, and for the presence of a dense 

£> . environment. 

ft 
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1. Introduction 

Snapshot fits of the broadband spectrum of gamma-ray burst (GRB) afterglows to 
standard forward shock models have been found, in many cases, to be consistent with an 
external environment density n < 1 cm~ 3 typical of a dilute interstellar medium (ISM) which, 
to first order, can be taken to be approximately independent of distance from the burst (Frail 
et al 2001). In other cases, the forward shock is better fitted with an external density which 
depends on distance as p oc R~ 2 , typical of a stellar wind environment (Chevalier & Li 1999, 
2000; Li & Chevalier 2003). The two types of fits have been critically analyzed by e.g., 
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Panaitescu and Kumar (2002), the conclusion being that at least some bursts may occur in 
high mass-loss winds, as expected from massive progenitors. The parameters for such wind 
fits are uncertain, due to poorly known stellar mass loss rates. 

In this Letter we show that observations of prompt optical/IR and/or X-ray emission 
attributable to reverse shock emission could constrain the GRB environment. The reverse 
shock emission tends to be in the regime where the electron cooling time is shorter than the 
dynamical expansion time. In high density environments, such as a stellar progenitor wind 
or a dense ISM in early galaxies, the self-absorption (SA) frequency is much higher than in 
the normal ISM, and it could be higher than the typical injection peak frequency ( Wu et al. 
2003). Here we argue, from general radiative transfer considerations, that in such situations 
when the emission is in the fast cooling regime and the SA frequency is higher than the 
injection frequency, the SA frequency and its scaling are different from, and the flux at the 
SA frequency is appreciably larger than, what had been previously estimated. This implies 
a different light curve time behavior for the afterglow prompt flash in a dense environment. 
This is of significant interest, since observations of the SA frequency and the net flux could 
provide constraints on the otherwise poorly known wind mass-loss rates of progenitors stars, 
or on the presence of a dense ISM. These new features are expected to be pronounced if the 
inverse Compton process does not play a dominant role. This could happen in the case of 
highly magnetized fireball ejecta, whose presence has been suggested by some recent studies 
(Zhang, Kobayashi & Meszaros 2003; Kumar & Panaitescu 2003) and by reports of high 
gamma-ray polarization (Coburn & Boggs 2003). The strength of the SA features would 
then provide a constraint on the magnetization parameter. 



2. The Model 

We consider a relativistic shell (fireball ejecta) with an isotropic energy E, an initial 
Lorentz factor r] and an initial width A expanding into a surrounding medium with a 
density distribution p = AR~ 2 . The shell width A is related to the intrinsic duration of 
the GRB as A ~ (l + z^cT (Kobayashi, Piran & Sari 1997) where z is the redshift of the 
burst. The interaction between the shell and the wind is described by two shocks: a forward 
shock propagating into the wind and a reverse shock propagating into the shell. The shocks 
accelerate electrons in the shell and in the wind material, and the electrons emit photons 
via synchrotron process. 

The evolution of the reverse shocks can be classified into two cases depending on the 
value of the initial bulk Lorentz factor 77 relative to a critical value r] cr = ((l+z)E/4:7rAc 3 T) 1 ^ 
(Sari & Piran 1995; Kobayashi & Zhang 2003b). If 77 > i] cr (thick-shell case), the re- 
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verse shock starts out relativistic and it significantly decelerates the shell material, to 
r ~ i] cr . If i] < i] cr (thin-shell case), the reverse shock is initially Newtonian and be- 
comes only mildly relativistic traversing the shell. Most bursts in a wind environment fall 
in the thick-shell class, (rj > i] cr ), because the high density implies a critical Lorentz factor 
f] cr ~ 70 C^^E^ 4 A^^T^^ lower than the typical value in fireball models rj > 100 (e.g., 
Lithwick & Sari 2001), where C = (l + z)/2, E 53 = E/10 53 ergs, the duration T 50 = T/50 sec, 
and A n _ 7 = A/5 x 10 11 g cm -1 is a typical wind mass loss rate. We are especially interested 
in long bursts with t x > 1 minute, for the purpose of analyzing the light curve behavior 
of the reverse shock emission. We focus on the thick-shell case, for which one can take the 
shock crossing time t x ~ T and the bulk Lorentz factor during the shock crossing (t < t x ) 
is T ~ r] cr . 

If we neglect synchrotron self-absorption (which is discussed in the next section) the 
spectrum of the reverse shock is described by a broken power law with a peak F u , max and two 
break frequencies: a typical frequency u m and a cooling frequency v c (Sari, Piran & Narayan 
1998). Assuming that constant fractions (e e and 6b) of the internal energy produced by the 
shock go into the electrons and the magnetic field, the reverse shock spectrum at a shock 
crossing time t x ~ T is characterized by (Kobayashi & Zhang 2003b), 

u e (T) ~ 1.5 x 10 11 (^^4lJ 2 Hz (1) 
u m (T) ~ 5.0 x 10 12 C- 1/2 <- 2 4 /2 -i^53 1/2 ^ii.7% 2 T 50 1/2 Hz (2) 
F v , max {T) ~ 95 d-^h^E^Ali^Tw 1 Jy (3) 

where e B -i = e^/10 -1 , e e _ 2 = e e /10~ 2 , i] 2 = 77/IOO, d = d L (z)/(2 x 10 28 cm), d L (z) is the 
luminosity distance of the burst, and g^(1) ~ 2 x 10 28 cm for the standard cosmological 
parameters (Q m = 0.3, Q\ = 0.7 and h = 0.7). The value of es(e e ) assumed for the reverse 
shock emission are larger (smaller) than in some previous broadband fits of afterglow forward 
shocks (e.g. Wijers & Galama 1999; Panaitescu & Kumar 2001). However, the values of 
the equipartition parameters can in principle differ in the forward and reverse shock regions, 
and as shown recently, 6b in the reverse (fireball ejecta) shock could be much higher than in 
the forward shock (Zhang, Kobayashi & Meszaros 2003; Kumar & Panaitescu 2003). 

3. Self-absorption and O/IR reverse flash 

In a stellar wind the external density at the initial interaction with the shell is much 
larger than in the typical ISM, hence the cooling frequency v c is lower than the typical peak 
frequency u m , and the synchrotron self-absorption (SA) frequency u a can be much higher 
than u m : v c (T) < v m (T) < v a (T). As suggested by some recent work (Zhang, Kobayashi 
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& Meszaros 2003; Kumar & Panaitescu 2003; Coburn & Boggs 2003), the fireball ejecta 
could be highly magnetized. In such cases the condition e e /e B <C 1 can lead to a Compton 
parameter F < 1 and the inverse Compton process is not important for electron cooling 
(the opposite case Y 3> 1 is discussed below). A consequence of the synchrotron dominance 
is that self-absorption suppresses the emission below the SA frequency u a , and prevents the 
electrons from cooling down to the Lorentz factor 7 C corresponding to the cooling frequency 
u c . The suppressed radiation energy is partly redistributed among the electrons, and results 
in a distinctive hump in the spectrum of the reverse shock emission. 

The reverse shock injects electrons with a power law energy distribution N^d'y oc 
r )~ p d r y{ r ) > 7 m ). The energy deposited in electrons with Lorentz factors between 7 m and 
7 a is E e ~ {p — l)N t Y^ m m e c 2 / {p — 2) where N t is the number of electrons in the shell, 
7 a is the Lorentz factor corresponding to the SA frequency, and we assumed p > 2 and 
7a 3> 7 m (so E e is essentially the total electron energy). This energy is redistributed among 
the electrons and photons in the optically thick regime on a timescale comparable to the 
cooling time (Ghisellini & Svensson 1989). Since the reverse shock is in the fast cooling 
regime (7 C < 7 m < 7 a ), the electrons have enough time to redistribute the energy. In a 
dynamical time, the energy E e is radiated as photons around u a . Assuming p = 3, the flux 
at u a is given by 

A-ndl v a T " v a 

A simple estimate of the maximal self-absorption flux is given by a black body flux with 
the reverse shock temperature (e.g. Sari & Piran 1999; Chevalier & Li 2000), F bb = 2tt(1 + 
z) 3 islm e T'y a (R±/dL) 2 where R± ~ 2rcT/(l + z) is the observed size of the shell. Equating 
Fj£ ~ F Ua , we obtain the S A frequency 

u a (T) ~ 8.1 x lO-C-3/M^/r 2£ Vu A 2/r E vi4 T -n/i4 ^ (5) 

An alternative derivation of the SA frequency is obtained by requiring the electron syn- 
chrotron cooling rate and heating rate (through absorption) to be equal at 7 a . The cross 
section for the synchrotron absorption process is approximately a s ~ 7 _5 r e ri (e.g. Ghisellini 
& Svensson 1991) where r e = q 2 /m e c 2 and r L = m e c 2 /q e B are the classical electron radius 
and the Larmor radius, q e is electron charge, B = (327te B T 2 pc 2 ) 1 / 2 is comoving magnetic 
field strength behind the shock. Using this cross section and the photon flux determined 
by eq. (4), we can evaluate the heating rate, while the cooling rate is given by the electron 
synchrotron power. Equating these rates reproduces the SA frequency of eq. (5). 

Following Kobayashi & Zhang (2003b), we obtain the scalings at t < t x of the spectral 
quantities v c oc t, v m oc t" 1 , v a oc t -5 / 7 , F v ^ max oc t°. Although v c itself is not observed, as a 
result of the self-absorption, the flux at v > v a can be still estimated from these scalings, 



F (T) ~ —— — - ~ fil£lH^ F (A) 
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because the electron distribution N e producing some observed frequency v (say in the K-band 
v k > v a) (and frequencies above this) is determined by the distribution of injected electrons 
at the shock, and by the synchrotron radiation cooling. Therefore, even though the flux at v a 
has a hump which was previously unaccounted for, we can apply the conventional synchrotron 
model to estimate the light curve at a frequency v > v a . The optical/IR luminosity initially 
increases as ~ {y jv a ) 2 F Va oc t 15 ^ 7 . When the SA frequency passes through the observation 
band v ohs at t a , the flux reaches a peak of F Uoba (t a ) ~ 2(u m u c ) 1/2 F Utmax /u obs and then it 
rapidly decreases. For electrons in the hump which are quasi-thermally distributed, the 
flux beyond the peak would decrease cx exp(— t 5 ^ 7 ) (or, if the emission of the quasi-thermal 
electrons above the peak is fitted by a power law v~ k above the SA frequency, the decrease 
is oc HsAXfc-i)^ the flux dropping by a factor 71 ~ 2[u obs /u m {t a )] ( - p - 2 ^ 2 . The time t a of 
the absorption peak passage through the K-band, the peak flux F UK (t a ) and the peak flux 
contrast 7Z relative to the subsequent power-law decay value are given by 

t ~ 20 r^ 10 e 2/5 e 1/10 A 2/5 rVio^-i/io -7/5 (&) 

F VK {t a ) ~ 1 d-'Cee^E^T^u-]^ Jy, (7) 

7? _ 7 .1/10-4/5 -1/5,-3/10^3/10 l T -3/10 -1/5 (R) 

~ ' S e e -2 e B ,-l A 11.7 -^53 % J 50 U obs,U-2 

where v ohs ^.i = u obs /(1.6 x 10 14 Hz) and p = 3 was assumed. At this turnover, a color 
change from blue to red is expected (see fig 1). Since the polarization is zero for an optically 
thick quasi-thermal spectrum, the reverse shock shell can emit polarized photons only above 
the turnover. If a turnover characterized by t a ,F VK (t a ) is observed, we can constrain the 
mass loss rate A, assuming that the redshift z is measured and the GRB explosion energy 
E is determined, e.g. from late time bolometric afterglow observations. The peak flux 
determines the equipartition parameter e e via eq. (7), and the peak time gives through eq. 
(6) a constraint on the mass loss rate, 

A ~ 5 x 10 11 eB ^C 3/4 e -i 2 ^3 1/4 T 5 f^ 2 i4,(^/20sec) 5 / 2 g cm - ( 9 ) 

where the dependence on the parameter t B is weak. 

After the optical/IR reverse shock emission drops to the level expected from the usual 
synchrotron model, the flux decreases slowly as ~ {v m / Vc^^iv / v m )~ p l 2 F v ^ max ~ t^(p~ 2 )/ 2 . 
Beyond a timescale comparable to the burst duration T, the optical/IR emission fades 
rapidly, because no further electrons are shocked in the shell (allowing the initially weaker 
but longer lasting forward shock to gradually dominate). The angular time delay effect 
prevents the abrupt disappearance of the reverse component, whose flux decreases steeply 
as ~ t^(p+ 4 )/ 2 (Kumar & Panaitescu 2000; Kobayashi & Zhang 2003b). No color change is 
expected around this break. The flux at this break for p = 3 is 

F VK (T) ~ 0.1 rf- 2 C 3/4 ^,-2i^i^l( 2 7^r^T 5 - 5 /V-^ 4 , Jy. (10) 
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These two breaks are schematically shown in Fig. 2, the break sharpness being an idealiza- 
tion; in reality these would be rounded. In our treatment, we have ignored pair formation 
ahead of the blast wave, e.g. Beloborodov (2002), which may modify the forward shock 
spectrum and light curve. 

4. IC effects and prompt reverse X-ray flash 

The above discussion assumed magnetized ejecta with e e /e£ <C 1 and Y < 1. However, 
for weaker ejecta magnetization with e e jts ^> 1, the inverse Compton (IC) process can af- 
fect the observed spectrum and the light curve, resulting in a reverse shock prompt X-ray 
flare. (For forward shocks, the importance of IC emission has been discussed, e.g. by Sari & 
Esin 2001; Panaitescu & Kumar 2001). Although the fraction of photons scattered to higher 
energies is small, they carry away a majority of the electron energy. The energy available 
for the synchrotron process is reduced from the injected energy E e by a factor of (1 + Y) 
where Y ~ (cg/e^) 1 / 2 is the Compton Y parameter (e.g. Zhang & Meszaros 2001). The IC 
process enhances the electron cooling, so v c is smaller by a factor of (1 + Y) 2 than its previous 
(synchrotron only) value. The thermal bump is shifted to a lower frequency, becoming less 
prominent or even disappearing. As a consequence of the IC cooling, the flux at v a becomes 
smaller by a factor of ~ (1 + Y) than the value given by eq. (4). If v a is shifted below u m , 
most of the energy available for the synchrotron process is radiated between v a and v m . The 
flux at v a is reduced from eq. (4) by a factor of ~ (1 + Y)( n f m / n f a ) ~ (1 + Y)(y m / V a ) 1//2 . By 
equating the flux at v a and the black body flux at the reverse shock characteristic temper- 
ature 7 a oc z/^ 2 , we can obtain the SA frequency. For the same parameters as for eq. (4) 
but taking e e _i = Ce/lO" 1 and e B _ 3 = e B /10" 3 , we get u a (T) ~ 5.8 x 10 13 C~ 3/U ^fli^-3 
A[ 7 ,EltT^ l,1A Hz and u m (T) ~ 5.0 x 10 13 C^-i^U^'* Mi.*^'* Hz ' To P ro " 
duce a significant bump, the SA frequency u a should be much higher than the typical fre- 
quency v m . When v a < u m , as in this case, the contrast 1Z is < 2 and the bump practically 
disappears. The optical/IR light curve initially increases as ~ {y jv a ) h l 2 F l , a oc t 5 / 2 , and after 
v a crosses the observation band, it decreases as t~( p ~ 2 )/ 2 . The transition is expected to be 
gradual and smooth. 

Since at the shock crossing time the forward and reverse shocked regions have roughly 
comparable energy, the vF v peaks of the synchrotron emissions reach roughly similar levels. 
Assuming that the characteristic reverse shock IC frequency v a ,ic ~ l 2 a v a is close to the 
typical (peak) frequency of the forward shock synchrotron emission, we can infer that in the 
case of Y 1 the reverse shock IC component is generally masked by the forward shock 
synchrotron emission, whereas in the case of 7 > 1 the reverse shock IC peak sticks out 
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above the forward shock synchrotron peak. Therefore, for weakly magnetized fireballs (with 
Y > 1), a prompt X-ray flare is expected from the reverse shock. The characteristic photon 
energy and the flux at this frequency are hu a jc(T) ~ 20 (~ 3 ^ 7 e^^e^ 1 ^ E^ A^jT^^ 7 

keV and F UaJC (T) ~ 1 rf-^ 31 / 1 ^^!^! 7 ^^ 14 ^!!^ 4 ^) 17714 m Jy> respectively. The typical 
duration of this X-ray flare at v > v a jc is of order the shock crossing time t x ~ T = 50T 50 
s. After the reverse shock crosses the shell, electrons are no longer heated and the (on-axis) 
synchrotron flux at v > u a {T) drops, as does the the IC emission at v > v a jci and one 
starts to observe high latitude emission. The X-ray emission from the reverse shock decays 
steeply as t~(p+ 4 )/ 2 , whereas the forward shock emission decays oc t~( 3 P~ 2 )/ 4 . Thus the slower 
decaying forward shock component eventually starts to dominate. 



5. Dense ISM in Early Galaxies 

We consider a specific model in which the ISM density of early galaxies scales with 
redshift as n ~ (1 + z) 4 n cm -3 (e.g. Ciardi & Loeb, 2000). In this case (or in general when 
the ISM density is much larger than a typical no ~ 1 cm -3 at z — 0), a discussion analogous 
to that of the previous section can also lead to a bump in the reverse shock spectrum and 
in the optical/IR light curve. 

The critical Lorentz factor classifying the evolution of the reverse shock is given by (e.g. 
Kobayashi 2000) r]' cr ~ 90 n ^Q^E^T'l'* where Ce = (l + z)/Q and T 2 . 2 = T/150 sec. A 
large fraction of GRBs are expected to be classified as thick shell cases, whose spectral quanti- 
ties are given by (Kobayashi & Zhang 2003a) v c (T) ~ 2.0 x l^n^^^e'^E'^T'l' 2 Hz, 

u m (T) ~ 7.2 x lO^Ce^K-i^ Hz and F,, max (T) - 27nl /4 D'Xi 1/4 ^ 2 _ 1 E%% 2 %:^ 
Jy, where D = d^(z)/(1.5 x 10 29 cm) is the normalized luminosity distance and D(z = 5) = 1. 
Equating Fj% ~ F Va , we obtain v a {T) - 3.7 x l^r^' 7 ^eJ^^B^^T^^ Hz. Using 
the scalings by Kobayashi (2000), one can show that u a oc £ -3 / 7 during the shock cross- 
ing. The optical/IR light curve initially increases as ~ (u/u a ) 2 F Ua oc t 9 ^ 7 . When u a passes 
through the K-band at t a ~ 5 e^ 2 e^_ x E\j^T 2 ^ 2 v~ll^ 2 sec, the flux reaches a peak 

of F UK (t a ) ~ Q.lD^ 2 C > Qe e ^2E^T2~2^ob s 14.2 Jy> and ^en it rapidly decreases by a (bump con- 
trast) factor of 1Z ~ 30 rig 1 ^ 4 C(T l72e ^-2 e B 1 -i 7 ?2 _lz/ o6s 2 i4 2 (P = 3)- After the emission drops to 
the level expected from the usual synchrotron model, the flux keeps a constant level, and 
then it rapidly fades as t~(p+ 4 )/ 2 after a time comparable to the burst duration. 
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6. Discussion and Conclusions 

We have analyzed the prompt afterglow emission from the reverse shocks of GRBs 
occurring in dense environments, such as the stellar wind of a massive progenitor, or a 
dense ISM as might be expected in early galaxies. Usually in the fast cooling case, the vF v 
flux is normalized at the typical frequency u m by using the energy ejected into electrons. 
However, here we point out that if the synchrotron self absorption frequency v a is higher 
than the typical frequency v m (and v c ), this usual prescription should be inapplicable for 
estimating the flux at and below v a . Such conditions can occur in the reverse shock of a 
highly magnetized fireball in a dense environment. The radiation flux suppressed by the 
self absorption effect is redistributed among the electrons and photons in the optically thick 
regime, and most of it is emitted at ~ v a in a dynamical time. As a result, the self-absorption 
frequency is different and scales differently with the shock parameters, and the flux at the 
self-absorption frequency shows a bump which is a factor 2(z/ a /z/ m )( p ~ 2 )/ 2 (~ several) above 
the usual power law flux estimate for typical parameters. The flux well above v a is the same 
as before, but the flux below v a is larger by the same factor 2(v a / v m )^ p ~ 2 ^ 2 . This results 
in a new type of temporal behavior for the prompt optical flash of afterglows from fireballs 
in dense (e.g. wind or early galaxy) environments. These new features will be prominent 
when the inverse Compton process is not important for electron cooling, i.e. for ejecta with 
Ce/e_B *C 1. This may be the case in fireball ejecta which are highly magnetized, as suggested 
by some recent studies. If, on the other hand, the burst occurs in a dense environment and 
such features are absent, this may be an indication that e e /e£ ^> 1, and in this case a prompt 
reverse shock X-ray flare is expected, which for a brief time dominates the forward shock 
but decays faster than it. 

Massive stars appear implicated in producing long gamma-ray bursts (as seen from 
the detection of a supernova associated with GRB 030329, e.g., Stanek et al 2003). Wind 
mass loss is expected from such stars previous to the GRB explosion, but snapshot fits to 
forward shock late emission (e.g., Panaitescu & Kumar, 2002) are compatible with such wind 
mass loss in only a handful of cases. In general the parameters of stellar winds are poorly 
known, and the uncertainties are further increased at high redshifts, where massive stars 
are expected to be metal poor. For this reason, signatures of a wind mass loss or a dense 
environment would be extremely valuable, both for GRB astrophysics and for tracing the 
properties of star formation at high redshifts. The prompt optical flashes expected after tens 
of seconds from the reverse shock in a dense environment would give characteristic signatures 
in the spectral and temporal behavior. These may help to test for the presence of winds 
and constrain the wind mass loss, in moderate redshift environment, or alternatively, at high 
redshifts they may provide evidence for a denser ISM than at low redshifts. In such winds or 
dense ISM the spectra and light curve time behavior can also give constraints on the strength 
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of the magnetic field in the ejecta. Large numbers of prompt X-ray detections with future 
missions such as Swift, complemented by ground-based follow-ups, should be able to test for 
such wind or dense ISM signatures and trace any changes with redshifts, if they exist, thus 
constraining the GRB environment as well as the radiation mechanisms. 
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Fig. 1. — Reverse shock spectrum in a dense environment or wind when synchrotron losses 
dominate: with self-absorption (thick solid) and without self-absorption (thin dashed). The 
schematic self-absorption maximum would appear as a rounded thermal peak. The previous 
self-absorbed flux estimate is shown by hashed lines. The correction factor 2(u a /u c ) ( - p ~ 2 ^ 2 oc 
£(p- 2 )/7 i s slightly larger at later times (t ~ t x ), the value of ~ 7 is evaluated at t a for typical 
parameters. 
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Fig. 2. — Schematic optical light curve for a synchrotron dominated fireball in a dense 
environment or wind: reverse shock emission (solid) and forward shock emission (dashed). 
The hashed line shows a previous estimate. Time scales are rough estimates for the typical 
parameters. 



